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ABSTRACT 



This report presents the results of an experimental 
investigation of the behavior of the turbulent early wake 
behind three flat plates of varying trailing edge config- 
uration submerged in water. 

The circumstances surrounding "singing" propellers and 
the influence of vortex formation upon this and other 
hydroelastic vibration phenomena are discussed. The sig- 
nificance of the vortex as a forcing function is emphasized. 

The influence of plate vibration upon the early wake is 
essentially eliminated in the investigation by rigidly 
clamping the test plates. 

A prediction of plate profile drag coefficient and of 
dimensionless shedding frequency, the Strouhal number, is 
made based upon von Karman's vortex street analysis. 

The longitudinal and transverse distribution of local 
mean velocity, turbulence intensity, and static pressure 
were measured in the early wake of three selected test 
plates of varying trailing edge configuration. Turbulence 
energy density spectra were obtained at selected positions 
in the early wake. Profile drag coefficients were computed 
and vortex shedding frequencies measured, and are compared 
with predicted values. 
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Ap = discrete area under vortex shedding frequency peak in 
turbulence energy density spectrum, dimensionless. 
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Op = drag coefficient. 
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C v = 1 . also vibration analyzer calibration constant, volts" 1 
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E(n) = spectral distribution of turbulence kinetic energy density 
per unit mass, ft 2 /sec 2 per cps. 
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X 



gage, feet of water. 

H-j-g = average effective total head seen by turbulence gage. 

H{.g = turbulent fluctuation of effective total head seen by 
turbulence gage. 

k = vortex strength, ft^/sec. 

E t g = turbulence gage sensitivity, mv/psi or mv/ft of water. 

K tp = transducer-preamplifier system sensitivity = GpK^.g. 
n = frequency, cps. 

An = band width, cps. 

n c = vibration analyzer band width center, or tuned, frequency. 

ng = Helmholtz resonant frequency. 

n Q = natural or resonant frequency. 

n v = vortex shedding frequency, cps. 

p = instantaneous pressure intensity, psl. 

p = average pressure intensity, psl. 

p' = turbulent fluctuation of pressure intensity. 

p 0 = free stream pressure intensity. 

Rd = vibration analyzer meter reading. 

Re = Reynolds number. 

Re c = chord length Reynolds number = Uc/v. 

Re^ = diameter Reynolds number = Ud/v. 

S = Strouhal number. 

S b = Strouhal number based on transverse vortex separation = 
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u = instantaneous local velocity in the downstream or 
x-direction. 
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I. INTRODUCTION 

A. General Statement of the Problem. 

One of the most vital problems in ship design today, and 
particularly in warship design, is the reduction of noise and 
vibration. Of the many and varied sources of noise and vibra- 
tion in ship structures, one of the more significant is the 
hydroelastic vibration of propeller blades and control sur- 
faces. The phenomenon of "singing" propellers has interested 
workers in several fields since the 1930's, when higher pro- 
peller speeds and the use of airfoil blade sections instead 
of ogival sections became prevalent ^ . Fundamentally, blade 
vibration of some sort has been held responsible in the 
literature for this phenomenon, but writers have differed in 
their analyses of the vibration and in their several explana- 
tions for the cause of this vibration^ 1- ® ^ . 

While propeller in-flow variations, blade back pressure 
variations, and wheel unbalance may indeed contribute to the 
phenomenon, most writers, however, agree that the probable 
cause of singing appears to be the formation and shedding of 
vortices in resonance with the natural frequency of the blade 
and its entrained water. If the vibrations set up by the pro- 
cess of forming of eddies at the trailing edge of a propeller 

blade are in the audible range of frequency, a singing note 

rql 

may be emitted 1 ”^ . It has been found that for a given pro- 
peller system singing occurs in a relatively narrow frequency 
band and over a relatively narrow range of propeller rpm. 

While audible singing is generally not encountered in 
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larger foil or plate structures such as rudders, hydrofoils, 
submarine diving planes and other control surfaces, vibration 
and flutter are significant problems ^ 11 ^ . It appears 
reasonable to consider that vortex induced vibration may be 
at least in part responsible for this. 

It has been found in practice that modifications to the 
geometry of propeller blade sections will alter the character- 

r 12 1 

istics of singing or eliminate it entirely. Den Hartog 1 J 

and others discuss the procedure of sharpening the leading 

edge of propeller blades and turbine blade runners to eliminate 

undesirable singing. Propeller sections with sharp leading 

edges have, however, been known to sing violently, and several 

authors ^ 1 ^ discuss the practice of blunting the trailing 

edges of blades to alleviate this condition. The periodic 

transverse forces due to alternate shedding of vortices incite 

a resonant mode of vibration in some part of the propeller 

r o l 

blade, and according to Saunders 1 J the remedy is to change 
the shape of the trailing edge in some manner so that the vor- 
tex formation is eliminated or that the vortices and circu- 
lation are stabilized. One simple method is to cut off the 
trailing edge square, so.that the separation points are "well 
defined and fixed in position;" another is to plane the 
trailing edge flat at some angle to the base chord. Various 
other trailing edge geometry modifications are discussed in 
the literature. 

Since practical experience shows that modification of 
the trailing edge will alter the response of a blade section, 
and since it is not reasonable that such modifications would 



materially alter the mass or the internal damping of the 
r 7 1 

system 1 ,J , it follows that it is the forcing function, i.e., 
the shedding of vortices, which has in some way been modified. 

A large body of literature exists on the subject of wake 
flows behind various types of bodies. The majority of it, 
however, deals with the wakes behind bodies of revolution and 
bluff bodies, and almost all of it treats the so-called fully 
developed wake well downstream of the body. In order to con- 
duct a study of trailing edge vortex phenomena at their 
origin, it becomes necessary to investigate the "early wake" 
or the region within about ten characteristic dimensions 
downstream of a foil or flat plate. Such a study should pro- 
perly consider the distribution of the energy present in the 
early wake and the variations of velocity, pressure, and turb- 
ulence intensity in this region. Furthermore, the plate bound- 
ary layer itself and the free shear layer where vorticity 
becomes pronounced may be expected to influence these phenomena. 
These considerations have prompted an investigation by the 
authors of the early wake behind three flat plates of varying 
trailing edge configuration at the variable flow water tunnel 
of the Massachusetts Institute of Technology Hydrodynamics 
Laboratory. 

B. Scope of the Investigation . 

As part of a general investigation of wake mechanics 
conducted at the Massachusetts Institute of Technology Hydro- 
dynamics Laboratory and sponsored by the David Taylor Model 
mi 

Basin, Huval 1 recently conducted an investigation of the 
early wake behind a fixed flat plate having a semi-cylindrical 



leading edge, a chord-to-thickness ratio of 8.0 and a 60° 
included angle trailing edge. This particular plate was one 
of a series studied by Ippen, Toebes, and Eagleson^ 1 ^ for 
the purpose of investigating hydroelast^ic behavior of flat 
plates with various trailing edge configurations. The authors 
chose three other plates from this systematic series for their 
investigation; their dimensions and geometries are listed 
below in Table Is 

TABLE I 

DIMENSIONS AND GEOMETRIES OF TEST PLATES 

PLATE LEADING EDGE THICKNESS CHORD TRAILING EDGE 
A semi-cylindrical l/V 2" ^ *0 0 = 90 

B semi-cylindrical 1/V' 2" £ ) 

C semi-cylindrical 1/V' 2" [ 

These three plates correspond respectively to Plates 12, 11, 
and IS of Reference £ li+] • 

The experimental program for this investigation provided 
for the measurement of velocity and pressure distributions 
and distribution of turbulence intensity across the wake at 
several longitudinal positions in the early wake of the test 
plates, and for measurement of the turbulence energy spectra 
at several longitudinal positions on the wake centerline and 
in the vicinity of the paths of shed vortices. 

Plates A and B were originally selected from consider- 
ation of vibrational amplitude responses as a function of 
restraint (Figure 29 of li|); these particular plates repre- 
sented respectively the lowest and one of the highest 
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response curves presented. During the course of the investi- 
gation, it became desirable to include an intermediate con- 
figuration, and Plate C was incorporated into the investigation. 
As in Huval's investigation, the authors decided to clamp 
the plate rigidly in the water tunnel test section in order 
to be able to study the wake energy with minimum contribution 
from plate vibrations. In order to facilitate future compari- 
son with Huval's results, it was decided to conduct measure- 
ments at a free stream velocity, U Q , in the neighborhood of 
10 feet per second. 
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II. REVIEW OF PREVIOUS WORK 



A. Vortex Excited Blade Vibration 

Of the large number of explanations postulated for 
singing propeller blades over the past 30 years, the sugges- 
tion that vortex formation and shedding might be the 
responsible mechanism appears to be one of the most plausible. 
Experimental work conducted in the United States and in the 
Soviet Union in recent years^®’ 1 ^ support strongly the con- 
tention that singing occurs when the vortex shedding fre- 
quency is in resonance with one of the natural modes of 
vibration of a given structure. 

Considering that at resonance singing frequency coin- 
cides with vortex shedding frequency, Gongwer^®^ presented an 
empirical modification of the Strouhal formulation, suggesting 
that the characteristic trailing edge dimension should include 
a portion of the boundary layer thickness. Introducing this 
factor, he obtained experimentally from his work with vanes 



and where Uj is the component of the free stream velocity 
normal to the trailing edge of the body. Krivtsov also 



number for a given trailing edge thickness. 

Several writers have been concerned with the fact that 




( 1 ) 



where b'= t + 




( 2 ) 



illustrates the essentially constant behavior of Strouhal 



the singing behavior of identical propellers may be quite 



different. However, when one considers the processes in- 
volved in propeller manufacture and balancing, it is not at 
all surprising that two propellers of the same class might 
have slightly different natural frequencies, and that the 

frequencies of the individual blades might vary somewhat, 
r 7 1 

Lankester 1 J points out that slight differences in natural 
frequencies between blades are not at all uncommon. This may 
increase the response or decrease the equivalent damping of 
one propeller with respect to another having blades of iden- 
tical natural frequency. If propeller singing is due to 
self-excited vibration, a small difference in response could 
make all the difference between singing and non-singing. 

Propeller singing has been shown to be dependent upon 
ITll 

water temperature and ambient pressure, indicating further 
the flow dependence of the phenomenon. 

B. Vortex Formation and Wake Periodicity 

Since Strouhal's work in the latter part of the 19th 
century, a large number of investigators have studied the 
vortex shedding of cylindrical and spherical bodies, and 
more recently of bluff bodies. While this work is of con- 
siderable interest in theory and in practice, and valuable as 
background, it is not immediately applicable to the case of 
thin flat plates and other blade sections, and will be treated 
only cursorily. 

The Reynolds number dependence of the Strouhal number, 

IT 16, 1 

suggested by Raleigh 1 J has been clearly demonstrated 
recently for circular cylinders by Kovasznay^ 1 * - *^ and by 
Roshko^®^. Roshko presents empirical expressions for 



( 3 ) 



for two Reynolds number regimes of flows 

S d = 0.212 (1 - 21.2/Re d ) Uo<Re d <l50 

and 

S d = 0.212 (1 - 12.7/Re d ) 300<Re d <5000 (ij.) 

For wakes behind various bluff bodies, Ippen et_ al. re- 
view the work of several investigators and indicate that 
the general range of Strouhal numbers based upon the largest 
characteristic dimension normal to the direction of the flow 
is from about 0.12 to 0.21, i.e., below the value of 0.2 
applicable to circular cylinders. 

Fage and Johansen in experiments on the early wake be- 
hind inclined flat plates showed that Strouhal number is 

dependent upon angle of incidence, a, increasing as a becomes 

r iq 1 

small 1 ^ . Expressing S in terms of the separation between 
free shear layers, S was found to be approximately constant 
and 2 ; 0.16. Both in Gongwer's experiments and in other ex- 
periments on oval struts it was found that the effect of 
small angles of incidence was to alter the amplitude of 
vibration without materially affecting the frequency. 

In none of the theoretical treatments of the problem 
is an explanation of the actual mechanism of vortex formation 
presented. The work of Betchov^^ suggests that a possible 
approach may lie in the treatment of boundary layer oscil- 
lations. 

Von Karman"s classical treatment of periodic wakes has 
crystallized and stimulated significant work In this field. 
His stability criterion for a vortex array in a wake has 
theoretical and practical significance and will subsequently 
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be utilized to predict shedding frequencies and to analyze 
free shear layer separation. This concept and Birkhoff's 
modification will be discussed further under "Theoretical 
Considerations." Von Kantian' s treatment makes no attempt to 
describe the mechanism of vortex formation and is in fact 
independent of the mechanism and of the geometry of the 
structure behind which it applies. 

C. Wake Structure 

The similarity principles of von Karman applied to the 
Reynolds equations of motion have proven useful in treating 
axially symmetric or two-dimensional turbulent wakes. Un- 
fortunately, similarity hypotheses apply to a state of 
statistical equilibrium of flow processes relatively far 
downstream in a wake, and cannot be considered applicable to 
the early wake. 

Roshko's experiments^ 1 ®^ and those of Cooper and Lutzky 
present results of velocity and turbulence distributions in 
the relatively early wake behind cylinders and bluff bodies 
respectively. 

r 22 1 

Birkhoff 1 J and others point out that contrary to the 
predictions of the Kirchoff wake model, actual large wake 
underpressures exist. 

D. Trailing Edge Geometry 

Ippen et jQ. point out the effect of trailing edge 

geometry on drag coefficients for various shapes. Schlich- 

r 1 

ting 1- - JJ shows that the drag coefficient can be related to 
the vortex street through momentum analysis, but notes that 
this procedure is applicable only relatively far downstream 
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where the static pressure at the measuring station is equal 
to the free stream pressure. The correction methods of Betz 
and of Jones ^ 2 ^ may be employed in the early wake region, 
and have yielded valid results behind airfoil sections as 
close to the trailing edge as 0.05 chord lengths. 

E. Turbulence 

The complex nature of wake turbulence and turbulent 
(or even periodic) wake energy precludes any detailed dis- 
cussion here. Reference [28] in particular and portions of 

references [29] through [31] are of interest for this inves- 

M-ai 

tigation as are the results of Huval's experiments . 

Ambient turbulence background measurements were made 

r-221 

for the water tunnel test section by Perkins J and by 
fl^ 1 

Huval 1 indicating relative mean turbulence levels of 
about 0.03 and 0.065 respectively. Assuming isotropic tur- 
bulence with a general power decay law, and ergodic process, 
Huval's measurements indicated that the free stream turbulence 
was largely in its initial stage. The character of the free 
stream turbulence spectrum obtained by Huval indicates 
anisotropy in the low frequency end; the high frequency range, 
however, indicates isotropic statistical equilibrium. Huval 
notes that the high frequency spectrum appears to follow von 
Karman's suggested energy variation with n [!3J # This 
is in agreement with the Kolmogoroff spectrum law cited by 
Hinze^ 28 ^ . 

Hinze points out that in all cases where the mean vel- 
ocity shows a gradient turbulence will be anisotropic. It 
should then be expected that in regions of a wake where such 



gradients exist anisotropic turbulence will obtain. This con- 
dition is expected to exist in the early wake behind an object. 

Roshko^-' 1 '®^ notes that the velocity fluctuation in the 
wake of a vortex- shedding cylinder displays a predominant 
frequency which is the shedding frequency? he points out 
further that this predominant frequency will have random 
irregularities imposed upon it. 

Distributions of turbulence intensity, mean velocity, 
and pressure in a turbulent wake are discussed by Hinze. 
Theoretical treatments of these distributions are in general 
predicated upon similarity assumptions which are well known. 
Such assumptions, however, are not expected to be Justified 
in the early wake behind a structure. Cooper and Lutzky^" 1 ^ 
found that certain predictions made on the basis of similarity 
hypotheses were not upheld by experimental evidence in a 
region relatively close behind bluff bodies. They note that 
this discrepancy was not resolved. 



III. THEORETICAL CONSIDERATIONS 



A. Analysis of the Vortex Street . 

From consideration of complex potential theory, it can 
be shown^^ that the complex potential function due to an 
infinite row of equally spaced vortices is given by 

w = H ln sln f 

from which it can further be shown that for two such rows 
parallel to each other and separated by a distance b, 

w - [in sin ^ (z - - In sin (z - ^ + ip-)j (6) 

where s 

a = longitudinal vortex spacing 
b = transverse vortex spacing 
k = vortex strength 

dW 

Since u - iv = , differentiation and separation of 

components yields an expression for vortex velocity in the 
longitudinal directions 

Uy = ^ tanh ~ (7) 

From stability considerations, von Karman showed^-^’^) 
that the street can be stable only when vortices in one row 
are opposite the interval center between vortices in the other 
row, and then only when the ratio b/a = 0.281. 

Rosenhead (cited in 25 ) shows that rigid boundaries 
equidistant from the mean centerline of the vortex street 
may affect the stability ratio. If h = the distance between 
the rigid boundaries, then as the ratio a/h varies from 0 to 
0.8l5» the stability ratio has a unique but continuously 
changing value, going from 0.281 to 0.256. For values of a/h > 



0.8l£ stability exists for a certain range of values of b/a, 
and for values of a/h > I.I4.19, the configuration is stable 
for all values of b/a. Birkhoff (cited in li+.) predicts a 
stability ratio of b/a = 0.35> from inertia considerations. 

The test section of the water tunnel used in this present 
investigation is rectangular, measuring 9 inches high by 1 - 1/2 
inches wide. The test plate is mounted vertically in the 
section. The dimension h between rigid boundaries in this 
case is then 7 - 1 / 2 ". Assuming that at the trailing edge 
b = plate thickness, t, and proceeding from the von Karman 
criterion a/b = 0.281 as a starting point, the ratio a/h = 
0 . 118 £. Following a reiterative procedure and assuming 
linear decrease of stability ratio with increase in a/h, the 
Rosenhead treatment indicates that the stability ratio for 
vortex streets in the current investigation would be b/a = 
0 . 271 . In actuality the transverse spacing is expected to 
be greater than the plate thickness due to wake spreading. 

The assumption of linear decrease in b/a is considered 
adequate in the region of a/h concerned. It appears, then, 
that a stability ratio for b/a may lie in the region from 
about 0.27 to about 0.28 for this investigation. 

After the manner of Xppen^" 1 ^ and Birkhoff it can 
be shown that considerations of vorticity diffusion from the 
boundary layer can lead to an expression for the strength, k, 
of the vortices in the wake. If be the rate of shedding 
vorticity from each side of an obstacle, then as an approx- 
imations 

K i = J £ u idy = J If u i d y = ^i 2 



( 8 ) 



where u ± is the velocity outside the boundary layer near the 
separation point. If fi = H/K^, where K is von Harman* s rate 
of vorticity discharge from each side of an obstacle, then /3 
is a measure of the anihilation of vorticity at the trailing 
edge. Both Prandtl and Roshko( cited in llj. and 27) have found 
experimentally that /3 = 0.5 for periodic wakes. Applying this 
value and assuming that ui = U, we obtain; 

K = 0.25 U 2 (9) 

Introducing this into von Harman* s expression for vortex 
strength, we obtain: 

1c _ Ka _ 0.25 aU /, n i 

k " tf - u v " ! - Uv/lf 
Introducing this into Eq. (7) we obtain: 

, Uy. 



jp (1 - jp) = 0.125 tanh trb/a 



( 11 ) 



from which: 



■ 0 - — 0 . 5 (l - \/l “ 0,5 tanh trb/a ) ( 12 ) 

Taking as the stability ratio b/a = 0.271 and operating 
on Eqs. (10) and ( 11 ) we obtain: 

u v /U = 0.095S (13) 

k = 0.277 aU 

From von Harman's relationship between profile drag, D p , 
and wake vorticity, where Dp = kp^- (U- 2 u v ) + and from 

Eqs. (13) we obtain: 

D p = 0.0728 paU 2 (li+) 

from which: 

°Dp = °- U 56 a/t 



(IS) 
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and reintroducing the stability ratio b/a = 0.271: 

% = 0-537 b/t (16) 

It should be pointed out that Ippen et. al . assumed a 
value of K/u 2 » O.lj. which is in consonance with empirical 
results cited by Birkhoff ^7}, if, however, the experimental 
result of j3 = 0.5 is accepted and Eq. (8) is considered valid, 
this would lead to a value of u-l/U = 1.265* 

r 2 V 1 

Birkhoff 1 ,J points out that it is a corollary con- 
sequence of the theorem of conservation of momentum that in 
viscous wakes vortex transverse spacing tends to increase; 
this is experimentally observed. For slightly viscous fluids 
it is further noted that transverse vortex spacing is essen- 
tially equal to the wake thickness a little behind the ob- 
stacle. Birkhoff notes that for circular cylinders the 
critical diameter Reynolds number above which the wake first 
becomes periodic is lowered by external turbulence. For flat 
plates parallel to the flow the critical wake-thickness Rey- 
nolds number is about 3 .l|J4-\/ Re l where Re^ = ; this is the 

same order of magnitude as in the case of a circular cylinder. 

Following Ippen’ s procedure, the plate thickness Strouhal 
number for vortex shedding can be estimated from Eqs. ( 13 ) 
and (15): 

(i - £><»/•> - 2^2 

Introducing Eq. (16) we obtain; 

St = 0.24s I 



(17) 
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from which: 




( 18 ) 



If wake thickness and transverse vortex separation are taken 
as being approximately equal to plate thickness in the early 
wake, then for t = 0 . 25 " we obtain a predicted vortex shed- 
ding frequency of: 



It should be pointed out that none of the foregoing 
takes into consideration trailing edge geometry. Furthermore, 
these developments do not consider any Reynolds number de- 
pendence upon either Strouhal number formulation or periodic 
wake structure. 



= 11.76 U (cps) 



(19) 
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IV. DESCRIPTION OF EQUIPMENT AND INSTRUMENTATION 



A. The Variable Flow Water Tunnel . 

A detailed description of the variable flow water tunnel 
is given in Reference [ 33 3 • Photographs of the tunnel test 
section are presented as Figs. Ia. and lb.; Fig. Ic. is a 
schematic diagram of the tunnel and its associated flow cir- 
cuit. The hydraulic characteristics of the water tunnel are 
presented in Appendix A of Reference [llj.]. Tunnel operation 
is possible with any one or a combination of three centri- 
fugal pumps. 

Considerations of minimum mean level of background 
turbulence and consistency of pump operation led the authors 
to conduct all measurements with Pump III providing the flow. 
Referring to Fig. Ic., the flow was routed from Pump III to 
the right through the 18" distribution circuit, back to the 
left past Pump I (which was isolated from the circuit), and 
up through the water tunnel inlet riser. It had initially 
been planned to operate with both Pumps II and III in oper- 
ation simultaneously, with the stop valve between them shut, 
but flow reversal through Pump II and problems in velocity 
control led to the decision to isolate Pump II from the 
circuit. 

From the expanding tunnel inlet riser the flow passes 
through four No. 16 mesh screens and a 33" long contraction 
into the test section. The plexiglass test section Is 36" 
long and rectangular in section, measuring 9" high and 7-1/2" 
wide. From the test section, the flow returns to the 
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FIGURE la. Variable Flow Water Tunnel 
Test Section, Showing Location of Amp- 
lifier and Preamplifier. 




FIGURE lb. Variable Flow Water Tunnel 
Test Section, Left Side. 
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CIRCULATION PUMPS and DISTRIBUTION CIRCUIT 

(PLAN) 




Pump I 

FIGUBE Io. Water Tunnel and Flow Circuit. 
(From Huval, "The Early Wake of a Fixed 
Flat Plate.") 
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storage tank. 

B. Test Plates . 

The test plates studied in this investigation are 
described briefly in Table I and in more detail in Reference 
[1I4.]. They were manufactured from 1/V aluminum stock and 
carefully machined. The plates spanned the 9" vertical 
dimension of the test section, and presented a 2" chord in 
the direction of the flow. The test plates were anodized to 
prevent corrosion and pitting! this process proved to be 
quite successful. 

C. Test Plate Mounting . 

The test plates were mounted using the spindles illus- 
trated in Figure 6 of Reference [1I4.]. The spindles were 
attached to the plates with set screws such that the center- 
line of the spindles was essentially in line with the 
center of radius of the semi-cylindrical plate leading edge. 

The spindles were clamped to prevent rotation. 

Each of the spindles passed through a roughly rectan- 
gular l" by 2-1/2" hole in the upper and lower test section 
walls. These holes opened further into larger cavities. The 
upper cavity was originally sealed with a solid cover platei 
because of the tendency of air bubbles to accumulate in this 
cavity, a vent was subsequently installed in the cover plate. 

The lower spindle passed through a gasketed circular hole in 
the bottom cover plate and was again braced below the test 
section. This arrangement is due to the earlier investigation 
reported in Reference [1I4.]. 
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D. Traversing Mechanism , 

In order to obtain data on transverse and longitudinal 

wake structures with the flow measurement devices to be 

described subsequently, it was necessary to provide for 

traversing in these directions. The traversing mechanism used 
Men 

by Huval 1 JJ was employed in this investigation. The in- 
strumentation strut, located just downstream of the test 
plate, passes through the test section side walls through 
slotted blocks which permit a longitudinal traverse. These 
are shown installed in Figs. Ila and lib. The instrumentation 
strut itself, illustrated in Fig. lie, can be traversed in 
the transverse direction through the side plates and blocks. 

Longitudinal positioning of the strut is accomplished 
by means of scales graduated in 1/6 1+ inch mounted on the 
slotted blocks. A transverse scale graduated in 0.002 inch 
permits transverse positioning. 

E. Flow Measurement Devices . 

The total and static head tubes and the turbulence 
pressure transducer employed in the wake surveys are described 
in detail in References [13] and [32]. These instruments 
permit measurement of the total and static pressures in the 
wake behind the test plate, and will be treated separately 
below. In addition to these devices, a contraction manometer 
provides for measurement of mean free stream velocity. 

1. Turbulence Pressure Transducer . 

A cylindrical disc-shaped barium titanate ceramic 
transducer, 0 . 12 $" in diameter and 0.100 1 ' thick, is mounted 
in a tapered water-filled probe tip illustrated in Fig. III. 




FIGURE Ila. Left Side View of 
Test Section, Showing Position 
of Test Plate and Transverse 
Strut Positioning Scale. 




FIGURE lib. Right Side View of 
Test Section, Showing Total and 
Static Head Leads and Transdu- 
cer Cable to Preamplifier. 





FIGURE lie. Photographs of Traversing System. 
(From Huval, "The Early Wake of a Fixed Flat 
Plate.") 
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The ceramic is located behind a thin silver-plated brass 
diaphragm,, The tapered probe tip ends in a 1/16" diameter 
with a 0.025" hole drilled in it to permit "point" pressure 
fluctuation measurements. Pressure variations are transmitted 
to the ceramic through a water column expanding from 0.025" 
diameter to approximately 0.221" diameter, and through the 
diaphragm. 

Considerable care must be taken to insure that no air 
bubbles are present in the water column. Experience has shown 
that the presence of even minute air bubbles will materially 
affect the response of the system by altering the response 
of the water column. To this end, deaerated water is used for 
filling the probe tip. Reasonable success was experienced 
by pumping deaerated water through the probe tip for several 
minutes with a syringe. 

It was found earlier that the transducer did not respond 
to static pressure changes ^ . Recent improvements, however, 
have improved the response of the transducer and made quasi- 
static calibration possible. This will be discussed further 
in Appendix A. 

2. Transducer Response . 

The instantaneous response of a total head tube aligned 
parallel to the direction of mean flow, in Huval's notation 
is given by: 

H tg =P+P' +(u+u ') 2 

“ Y 2g 

= p + p' + u 2 + 2uu* + u’ 2 (20) 

Y Y 2g 2g 2g 

assuming that the terms involving v’ and w" do not significantly 



contribute to the total instantaneous head measured in the 
longitudinal direction. 

Now let e tg = the instantaneous output signal, of the 
turbulence gage in millivolts, and 
K-kg = turbulence gage sensitivity in millivolts 
per foot of water. 

Then the turbulence gage output signal becomes: 

®tg = K tg H tg = K t-g + P' + u 2 + uu‘ + u ,gN \ (21 ) 

\Y Y 2g g 2g / 

Now let Gp = preamplifier gain, and 

ep = voltage output of the preamplifier in milli- 
volts, 



Then : 

e p = GpE tg / p + p* + u 2 + uu' + u 1 (22) 

\Y 7 “ 2g “ 2g J 

Furthermore : 

e p = e p + e p ' 

*p = G P E tg ( p + 5: + u ' 2 ^ 

\ Y 2g 2g / 

Referring to the electronics equipment block diagram 
(Fig. V) it is noted that a DC filter (i.e., a 1 MFD capacitor) 
is installed between the preamplifier and the amplifier, since 
only the fluctuating components of the turbulence gage signal 
are of interest at the measuring instruments. Therefore, the 
input signal to the amplifier becomes: 

i 

«p' = «p “ e p 

= G d K + _ /p c + uu' + u'2 - 

P S V“ — — 2g ) (23) 

Letting e m = amplified turbulence signal input to the measur- 



ing instruments, and 
G a = amplifier gain, then 



(23) 



e m = G a G p E tg ( Pi + + u ' 2 ~ u>2 

\ Y g 2g 

Here e m represents the amplified turbulence gage signal in 
millivolts at the input to the measuring Instruments (rms 
voltmeters and vibration analyzer) . 

The true rms voltmeter output (measured by the thermo- 
couple voltmeter) becomes: 



) 



V rms = = G a G p E tg 



^p' + uu 1 + u' 2 - u' 2 j 2 
\Y ~ 2g / 



£ (24) 



Expansion of the terms within parentheses in Eq. (24) results 
in: 



f p ' + uu' + u'2 - u'^V s u 2 u' 2 + uu' (u 1 2 - u'2) 

Vy g 2g ' g 2 g2 

+ (u'2 _ u'2)2 + 2uu'p* 

4g2 Yg 

+ (u'2 - u® 2 )p ' + p*2 
Yg Y 2 

The full expression for rms voltmeter output becomes: 






rms 



e a s p K tg ( 

\ s 2 



® 2 + u u“3 + u 



• T-V « 



+ 2uu'p 
Yg 



u*2 p « 



Yg 




(25) 



Now if it may reasonably be assumed that u*«u, say for 

example that V\p7* ^C.0.1 u„ and that p'/Y is of the order 

of u 1 2 . only the first term within the parentheses makes a 
2g 

significant contribution to the rms signal. Hence: 



V 

v rms 



G a G p E tg ^^u'2 

s 



(26) 
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3. Total Head and Static Pressure Tubes . 

The total head tube employed is in fact a standard Pitot- 
static tube, but only the dynamic tap is utilized. Because of 
the longitudinal pressure gradients present in the wake of the 
test plate and their effect upon measurements, the static taps 
of the total head tube could not be used, and a separate static 
head tube was employed. The pressure taps of the static head 
tube were in line with the tips of the total head tube and the 
transducer probe. See Pig. III. Huval^*^ gives a detailed 
description of these instruments. Both pressure tubes are 
connected to mercury manometers. 

4. Contraction Manometer . 

Connected to the upstream and downstream ends of the water 
t unn el contraction are a mercury manometer and a manometer 
filled with blue-colored Meriam No. D-8325 indicating fluid 
(specific gravity = 1.75). These manometers and their piezo- 
meter taps are described in Appendix A of Reference [l4] . The 
contraction manometer was used to duplicate flow conditions in 
the water tunnel on separate days of operation. Calibration of 
this manometer was necessary for the particular flow conditions 
experienced, and is discussed subsequently. 

F. Electronics Equipment . 

Electronics equipment utilized in the measurements consis- 
ted of an electrometer preamplifier, a KinTel Model 11 IBP 
amplifier, a General Radio Type 762-B vibration analyzer, a 
Ballantlne vacuum tube voltmeter and a Sensitive Research Instr. 
Corp. Model A thermocouple voltmeter, and a Dumont Type 333 
dual-beam cathode ray oscillograph. A Hewlett-Packard Model 
200C audio oscillator was employed in equipment calibration. 



STATIC 

PRESSURE 

TUBE 





SECTION A-A 

ENLARGED SX 




SECTION B-B 

ENJWGED IOX 



TURBULENCE 

PROBE 



TOTAL 

HEAD 

TUBE 




FIGURE III. Details of Turbulence Probe and Pressure 
Tubes. (From Huval, "The Early Wake of a Fixed Flat 
Plate . " ) 
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These components are Illustrated in Pigs. IVa and IVb. 

(At the end of the investigation a General Radio Model 
1445-A vibration analyzer became available; this instrument 
extended the frequency range of spectrum measurements and 
was used in one location in the wake of Plate C for comparison 
purposes. ) 

It had originally been Intended to record the amplified 
signal from the pressure transducer on magnetic instrumentation 
tape for remote spectrum analysis on another wave analyzer. 
Recording difficulties made this impossible during the course 
of this investigation. 

A block diagram of the electronics components is pre- 
sented as Pig. V. 

1 . Electrometer Preamplifier. 

The primary purpose of the electrometer preamplifier is 
to provide high input impedance to the first stage of ampli- 
fication. Huval ^ describes the circuitry of this component, 
and points out that its addition to the system should permit 
response to frequencies below 1 cps. The measured gain of this 

i 

component (see Appendix A) is 0.4 over the majority of the 
operating frequency range. 

2. Amplifier . 

The EinTel amplifier used in the system has a very flat 
response over the operating range of frequencies. The design 
gain accuracy with a fixed gain setting is ± \% from DO to 2 ECS. 
The measured response characteristics (see Appendix A) over the 
operating frequency range point this out further. A DC filter 
is connected to the input of the amplifier. 



30 




FIGURE IVa. Electronics Equip- 
ment: Analysis Components. 




FIGURE IVb. Electronics Equip 
ment: Preamplifier and Ampli- 
fier Installation. 
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E 



FIGURE V Electronics Equipment Component Block Diagram 
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3. Wave Analyzer . 

The General Radio Type 762- B vibration analyzer is a 
constant percentage band width tunable voltmeter with a 
logarithmic scale. Its operating range is from 2.5 to 750 
cps. For the sharp selectivity position (in which all measure- 
ments were made) the relative attenuation is 3 db at a frequency 
differing by \% from tuned frequency. Calibration was per- 
formed using a signal from a Hewlett-Packard audio oscilla- 
tor, and is discussed further in Appendix A. A damping 
capacitance was inserted in parallel with the analyzer meter 
to reduce needle fluctuations. 

The Model 1554-A analyzer used briefly at the end of the 
investigation is a similar instrument, but has an extended 
high limit of 25,000 cps for its operating frequency range. 

4. Voltmeters. 

The Ballantine VTVM is a true rms-indicating meter only 
for sinusoidal waveforms. The thermocouple voltmeter, on the 
other hand, will indicate true rms voltages for any signal 
input . 

5. Oscilloscope . 

The Dumont dual-beam oscillograph was used primarily for 
visual monitoring of the amplified turbulence signal. In this 
respect it proved very useful in determining the presence of 
air bubbles in the turbulence probe tip. Such bubbles attenuate 
severely the higher frequencies of pressure fluctuation. See 
Figs. Via, VIb, and Vic. 
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FIGURE Via. Turbulence Signal 
Presentation in Wake of Plate 
0, y/t=0.5* Illustrating Pre- 
dominant Vortex Frequency. 



FIGURE VIb. Turbulence Signal 
Presentation on Wake Center- 
line, Illustrating Predomin- 
ant Double Shedding Frequency. 




FIGURE Vic. Turbulence Signal 
Presentation with Air Bubble 
in Transducer Probe Tip, Same 
Location as Figure VIb above. 
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V. CONDUCT OF THE INVESTIGATION 

A. Ambient Turbulence . 

A study was made of the ambient mean level of turbulence 

as indicated by the amplified rms voltage from the pressure 

transducer for various modes of pump operation. As noted 

earlier, Pump III and the previously described flow circuit 

were used in the final investigation; this combination 

yielded one of the lowest mean levels of turbulence. 

It had been decided to operate the tunnel at a flow 

velocity of approximately 10 fps to facilitate correlation 

ri?i 

with Huval's results 1 JJ . It subsequently developed after 
final calibration of the contraction manometer that measure- 
ments had been taken at a value of free stream velocity in 
the test section, U Q , of about 9*5 fps. 

Energy spectrum measurements were made for the empty 
test section at one longitudinal location on the tunnel 
centerline. 

B. Wake Survey Measurements. 

Pressure, mean local velocity, and turbulence intensity 
distributions across the wake were desired at representative 
locations in the wakes of the three test plates. A coordinate 
system was defined as follows; the x-axis was taken as the 
flow centerline, which was assumed to coincide with the plate 
and tunnel centerlines. The +x-direction was taken in the 
downstream direction, measured from an origin located in the 
transverse plane of the plate trailing edge. The origin of 
the y-axis was taken at the particular longitudinal station 
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where measurements were desired (x/t = 1.^0, 3-25, 6. 50. and 
9 . 75 ). The +y-direction was taken in the transverse (hori- 
zontal) direction toward the right side of the test section 
(see Fig. la.), it was decided to measure total and static 
head and turbulence rms voltage across the wake in increments 
of 0.1 plate thickness. Full wake surveys were made from 
y/t 1.5 to y/t - +1.5 in this manner at longitudinal 
stations x/t - 3.25 and 9 . 75 . Half -wake surveys (from the 

flow centerline to y/t = +1.5) wer e made at x/t = I.50 and 
6.50. 



C. Wake Energy Spe ctrum Measurements . 

Turbulence rms voltage and vibration analyzer measure- 
ments were taken at y/t = 0 (flow centerline) and y/t = +0.5 
(in line with the plate edge) at each of the longitudinal 
measuring stations mentioned above, for each of the three 
test plates. The empty test section background spectrum 
mentioned in Section V.A. above was taken at y/t = 0 and 
x/t = 3 . 25 . 



~ ^ngi tudinal Variation Meaaur^pn^ 

In addition, it was desired to study the longitudinal 
variation of pressure, velocity, and turbulence intensity 
along the centerline and at y/t = +1.5. T o this end the 
measurements described in Section V.B. were also conducted 
at these transverse locations at the following longitudinal 
locations s x/t = 0. 5 , 0.75. 1.0. 1.5. 2.0. 2.75, 3.25, 4.0. 
5*0, 6.5. 7.5, 8.5, and 9.75. 

E. Calibr ation of Equipment . 



Calibration of equipment is discussed in Appendix A. 
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F. Flow Duplication . 

With the foot valve at the discharge end of the water 
tunnel set and left in one position, it was found possible 
to adjust the pump discharge valve and the distribution 
circuit stop valve to duplicate very closely both velocity 
and pressure conditions within the test section. The Meriam 
fluid contraction manometer was used for this purpose. This 
procedure was followed in establishing flow during the 
measurements made for Plates B and Cj minor changes were 
made in the foot valve setting during the Plate A measure- 
ments. 
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VI. REDUCTION OF DATA 



A. Mean Flow Velocity, Up . 

The contraction manometer was used to measure free stream 
velocity. The Meriam blue fluid manometer lent itself to 
more precise reading and was thus used primarily rather than 
the mercury manometer. The contraction manometer was cali- 
brated against a Prandtl-pitot tube inserted through the upper 
tunnel wall at the upstream end of the test section. Velocity 
of the free stream was obtained from the expression: 



where Ah^ = Meriam blue fluid manometer difference 
reading in inches. 

B. Local Mean Velocity, u . 

Local mean velocity was computed from the total and 
static head tube manometer readings using the expression: 



Ah s = static head manometer difference in feet 

Az = difference between total and static pressure- 
side manometer readings, in feet 

Because of small variations in values of U Q for dif- 
ferent sets of data, it was decided to normalize u with respect 
to the value of u at y/t = +1.5 in every case. 

C. Turbulence Intensity, >y/u"»z 

Following the development set forth in Section IV, the 
turbulence intensity, u'2 , was calculated from the rms 



U 0 — 2 . 06/y/~SHp^ 



( 27 ) 




( 28 ) 



where: Ah^ = total head manometer difference in feet 



turbulence voltage measurements using the expression: 
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V u i 2 = V gns # 



(29) 



u GaGpEtg 



Turbulence intensity was normalized with respect to the 
value of u obtaining at the point of observation. 

D. Static Pressure, p. 

Static pressure was calculated from the static head 

tube manometer reading for each measurement location. It 

was normalized in terms of a pressure coefficient P ~ Po 

|PU 0 2 

where p 0 is the free stream pressure and U Q is the free 
stream velocity. 

E. Normalized Turbulence Energy Density, E(n)/ u' 2 . 

Normalized turbulence energy distribution density was 
calculated directly from the thermocouple voltmeter reading 
and the reading of the vibration analyzer at each observation 
frequency from the expression: 



This normalization process eliminates variations in 
absolute magnitude of energy density due to variations in 
u existing over the entire investigation. 

F» Profile Drag Coefficient, ^Dp. 

The profile drag coefficient was calculated at longi= 
tudinal stations x/t = 3.25 and x/t = 9.75, using the 
expression: ^ 



where the integration was performed numerically by the 
trapezoidal rule between the limits y/t = -1.5 and 
y/t = +1.5. 




Rd 2 (n) 



(30) 



n 




( 31 ) 



39 



Applicability of this expression to the early wake was 

r 2^1 

checked using the Jones Method J where? 







where: g 2 = total pressure at the measuring station 

P 2 = static pressure at the measuring station 
P 0Q = static free stream pressure 
q QO = total pressure minus static pressure in 
the free stream. 

Integration was performed numerically as above. 

It should be noted that the drag coefficient in terms 
of wetted surface is given by: 
c Dp = 16 CDwg 

for a chord-to-thickness ratio of 16 as is the case for 
these test plates^-^. 



VII. RESULTS OF THE INVESTIGATION 



A. Wake Structure. 



The results of the study of the wake structure of the 
three test plates are presented graphically as Appendix B 
to this report. The salient points evidenced by these plots 
will be discussed in this section. 

Figs. B-l through B-3 are summary plots, showing for 
each of the three test plates all of the quantities investi- 
gated in the wake surveys. Necessarily drawn to a reduced 
scale, these plots are intended merely to indicate in one 
location the general pattern of the normalized mean velocity 
profiles and longitudinal distributions, the transverse and 
longitudinal distributions of normalized turbulence intensity 
and the transverse and longitudinal distributions of normal- 
ized static pressure. 

Figs. B-4 through B-6 show for each individual plate 
the development of the dimensionless velocity profile. 

Values of u are normalized in each case with the faired 
value of u obtaining at station y/t = +1.5° Profiles are 
shown for positive values of y/t only. (This is the case for 
all remaining wake plots with the exception of the wake 
symmetry plots.) None of the plates show any distinct in- 
dication of rapid wake spreading. If the limits of the wake 
are arbitrarily supposed to be the value of y/t for which 
l { y/^| = ^ is seen kkat f° r Plate A the limits of the 



dl 
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wake gradually increase from about y/t = 1.0 to y/t = 1.2 
between the four longitudinal stations. For Plate B the wake 
spreads from about y/t - 0.8 to y/t = 1.2j while for Plate C 
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the wake spreads from about y/t = 1.0 to y/t = 1 .Ip. The 
velocity defects for Plates A and C at x/t = 3*25 are seen 
to be markedly greater than for Plate C, and the decay of the 
velocity defect Is more rapid for these two plates. The 
transverse distributions of velocity defect are essentially 
Gaussian. 

The increased velocity defect at the early stations for 
Plates A and C is illustrated further in Figs. B-7 through 
B-10, which present a comparison of the dimensionless velocity 
profiles at each of the longitudinal measuring stations. 

Close to the trailing edge both of these plates, which have 
relatively sharp separation points, have large defects, al- 
though at x/t = 9*75 all three profiles are essentially the 
same. Plate A, with its 90° reentrant angle trailing edge, 
is seen to be intermediate between the other two plates. 

Figs. B-ll through B-llj- present the transverse dis- 
tribution of the normalized turbulence intensity for all 
three plates at each of the longitudinal stations. The wake 
behind Plate C is seen to be the most turbulent in the 
centerline areaj that of Plate B is the least turbulent. 

Again Plate A is intermediate. Turbulence intensity decays 
rapidly downstream and at x/t = 9»7f> the intensity levels 
are essentially equal and approaching a constant value across 
the wake, indicating a decay to isotropic turbulence. Tur- 
bulence intensities appear to approach the free stream level 
by y/t = l.f> at each station for each plate, and the free 
stream level of turbulence intensity is seen gradually to 
increase as energy is diffused from the wake into the stream. 



Fig. B-15 through B-l8 illustrate the transverse dis- 
tributions of static pressure coefficient in the wake of the 
three plates at each of the longitudinal stations. Plates A 
and C with their more sharply defined separation points are 
seen to have extensive underpressure regions at x/t = 1.5 j 
P late B with its rounded trailing edge does not. As we move 
downstream the Plate B distribution becomes partially 
negative, and all three plots move together. At x/t = 9.75 
(see Fig. B-l8) the distribution plot for Plate A is seen to 
become markedly negative again. This particular presentation 
is not substantiated by the investigation of longitudinal 
pressure coefficient distribution illustrated in Figs. B-21 
and B-22, and may be due to data error or computation error. 
However, the possibility that this plot is correct must not 
be entirely discounted in view of the unusual trailing edge 
configuration of this plate. The Plate A plot in Fig. B-15> 
shows sharp excursions which may be not entirely due to 
scatter. That the sharp drop in pressure coefficient at 
y/t = 0.£ may not be due to data scatter is suggested by the 
fact that Plate B also indicates a smooth decrease in this 
region. It will be seen later that Plate A exhibits a very 
sharp longitudinal pressure coefficient gradient in this 
region and that reverse flow is also observed. 

Fig. B-19 presents longitudinal distributions of normal- 
ized mean local velocity at the centerline and at y/t = 1 .5. 
In these plots u is normalized with respect to free stream 
velocity, U Q . Reverse flow was observed for Plate A in the 
region from x/t - 0.75 to about x/t = 1.25» and a sharp 



longitudinal velocity gradient exists. Reverse flow was also 
observed for Plate C close to the trailing edge. 

Pig. B-20 shows longitudinal distributions of turbulence 
intensity for the three plates at the wake centerline and at 
y/t = l.f>. Again Plate C indicates the highest level of wake 
turbulence and Plate B the lowest; Plate A is intermediate. 

Pigs. B-21 and B-22 illustrate the longitudinal static 
pressure distributions for the wake centerline and at y/t = 
l.£ respectively. Fig. B-21 clearly shows the large under- 
pressure in the wake close to the trailing edge of Plate A. 
Plate C also demonstrates considerable early underpressure, 
while Plate B underpressure is relatively minor and very 
rapidly becomes an overpressure for a short period, as does 
the plot for Plate A, The centerline pressure for Plate C 
is seen to remain negative throughout the longitudinal range. 

Figs. B-23 through B-25 illustrate the essential 
symmetry of the wakes behind the three plates. Although some 
data scatter is seen, it appears that the exact centerline of 
the flow was slightly on the positive side of the established 
axis, fortuitously insuring that the flow centerline was 
actually observed during the half -wake surveys as well as in 
the full-wake surveys. 

B. Profile Drag Coefficients . 

Profile drag coefficients were calculated for the three 
test plates at longitudinal stations x/t = 3.2£ and x/t = 9.75* 
Calculations were initially made using Eq. (31 )» and then 
using the Jones method as indicated in Eq. (32). Result of 
these calculations are listed as follows in Table II. 



TABLE II 
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Profile Drag Coefficients 



PLATE 


x/ 1 


% [Eq. (31)] 


°D p [Eq. (32)] 


A 


3.25 


0.443 


0.141 


A 


9.75 


0.378 


0.171 


B 


3.25 


0.519 


0.351 


B 


9.75 


0.318 


0.195 


0 


3.25 


0.528 


0.196 


0 


9.75 


0.393 


0.308 



The above values of Cpp computed from Eq. (31 ) are 
corrected for the longitudinal pressure gradient by alge- 
braically subtracting from the momentum term a pressure 
coefficient defect term. This term was obtained by numer- 
ically integrating the defect of the pressure coefficient 
plot where the defect was taken as zero at y/t = +1.5. The 
results compare reasonably with the value of Cp p = 0.537 
predicted by Eq. ( 16 ). The results yielded by Eq. (32) are 
considered Incorrect. 

C. Vortex Shedding Frequency . 

Vortex shedding frequencies were measured by observing 
the peak responses in the normalized energy density spectra 
for the three test plates (see Appendix C). These shedding 
frequencies were observed most strongly at y/t = 0.5s in line 
with the edge of the test plate. The double shedding frequency 
is more clearly noticeable in the spectra taken at y/t = 0, 
on the centerline. In some cases the triple shedding frequency 
is also noticeable. 

Plate thickness Strouhal numbers were calculated for 
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each of the three plates, based upon averaged flow velocities 
and observed shedding frequencies; the results of these cal- 
culations are presented In Table III. 

In order to compare the results of these calculations, 

It was decided to normalize the Strouhal numbers with respect 
to Plate C. It was reasoned that for a small Reynolds number 
range the Strouhal number might be expected to remain 
essentially a constant. Since variations In Strouhal 
number must be due to variations in the transverse separation 
between vortices, Plate C was selected as a basis for com- 
parison since the square trailing edge provides distinct 
separation points in comparison with either of the other two 
test plates. It was then assumed for the sake of comparison 
that the transverse separation between vortices is indeed 
equal to plate thickness for Plate C. The Plate C Strouhal 
number was then taken as being the constant value for all 
three plates, neglecting small variations due to slightly 
different plate thickness Reynolds numbers, and this Strouhal 
number was considered a vortex-separation Strouhal number, 

Sb = nb/U. From this, the relative values of vortex separ- 
ation were calculated for Plates A and B. These results are 
also listed in Table III, and their physical significance 
is illustrated in Fig. VII. 



TABLE III 
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Strouhal Numbers and Relative Transverse Vortex Spacing 



PLATE 


o, _ nt 
s t " IT 


n 


u 


<3 _ nb 

S b “ IT 


b 


A 


0.237 


109.87 


9.66 


0.223 


0.236 


B 


0.251 


Hi].. 87 


9.53 


0.223 


0.222 


C 


0.223 


103.00 


9.63 


0.223 


0.250 




FIGURE VII. Relative Transverse Vortex Spacing 

The qualitative representation of vortex spacing in Fig. 
VII illustrates expected behavior in the wake of Plate B; 
the separation points move toward the centerline of the plate 
for a semi-cylindrical trailing edge. When it is realized 
that vortices do not appear at the separation points as dis- 
crete entities, but are gradually formed from vorticity 
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which diffuses out of the boundary layer, even the behavior 
of the vortex locus of Plate A appears plausible. A "cavity" 
of relatively large under-pressure exists immediately behind 
the reentrant wedge, as was pointed out in Section VII. A. and 
illustrated in Appendix B, and the flow would be expected to 
be drawn toward the centerline of the wake. The loci of 
vortices are illustrated as gradually diverging from the 
wake centerline in all three cases in conformance with the 
indication of wake spreading seen in Appendix B. 

Observed shedding frequencies were compared with the 
frequencies predicted from Eq. (19). This comparison is pre- 
sented below in Table IV. Again, exact agreement is not 
realized, but a reasonable approximation is indicated. 



TABLE IV 



Comparison of Actual and 
Predicted Shedding Frequencies 

P^ATE 



n. 



A 

B 



v actual 



1 10 



n 



v predicted 
Eq. (19) 
113.6 
112.0 



C 



115 

103 



113.1 
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D. Turbulence Energy. 

The results of the investigation of the turbulence 
energy contained at various stations in the early wake behind 
the three test plates are presented graphically in Appendix C. 

Since minor variations existed in the flow conditions 
obtaining for the individual locations, it was felt that no 
precise comparison could be made of absolute magnitudes of 
energy density. To this end, the spectra were normalized as 
indicated earlier. The data presented on these plots is 
independent of the local velocity at the point of measure- 
ment. 

The area under these spectrum curves is dimensionless 
and is a measure of the total energy present at a particular 
location in the wake. Integration of the energy density 
spectrum from zero to infinity should properly result in a 
value of unity; this is not the case for any of the spectra 
obtained. This indicates that a certain percentage of the 
total energy present is either not measured or not accounted 
fox in some manner. Possible explanations for this fact are: 

(1) The limits of integration are the finite limits 
of the wave analyzer. Clearly some energy must be present 
below 2.5 cps and above 750 cps. Subsequent measurement of 
one of the spectra with a General Radio Type 1554-A vibration 
analyzer indicates a considerable amount of energy present 
at about 1400 cps, which is possibly attributable to Helm- 
holtz resonance effects. (See Appendix D. ) It is not believed 
that any significant amount of energy is present below the 
low frequency limit of the wave analyzer because of the 
small frequency increment under consideration. 
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(2) The spectral measurements made with the wave anal- 
yzer are normalized with respect to measurements made with a 
true rms voltmeter. It is possible that the thermocouple 
voltmeter, being an all-pass energy measuring device and 
Indicating a true picture of the total energy present at a 
given location, measures components which the wave analyzer 
meter does not. The wave analyzer is a tunable voltmeter 
designed to measure Fourier components of vibration, and some 
energy components may not be recognized by it. During the 
course of the energy spectrum measurements, simultaneous 
readings were taken with the Ballantine VTVM. On the assumption 
that the foregoing might be correct, the energy spectra were 

i 

also normalized with respect to the Ballantine voltage, since 
it was reasoned that this instrument is similar to the meter 
installed in the wave analyzer. The results of these calculations, 
which are not plotted, still indicate an energy defect, but a 
much smaller and more reasonable one. 

Further investigation will be necessary to explain this 
discrepancy fully before any absolute magnitudes can be 
assigned to energy data. 

Because of the foregoing, no magnitude comparison of 
discrete energy represented by the vortex shedding frequencies 
could be made. In order to examine the vortex energy as a 
function of trailing edge configuration then, the ratio 
Ap/A g was studied, where Ap is the area under a discrete 
shedding frequency Cor double shedding frequency) peak and 
A g is the total area under the spectral plot obtained. This 
ratio is then a measure of the percentage of energy present 
in the vortex shedding phenomenon for each of the spectra. 
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The ratio A p /A g is greatest at any given location for 
Plate B and least for Plate 0; Plate A is again intermediate. 

This fact was noted for both the fundamental and double 
shedding frequencies. The ratio A p /A g shows a linear decay 
downstream for both Plates B and C. The value for Plate A 
increases to a maximum in the vicinity of x/t = 3.25 and 
then exhibits essentially a linear decay with x/t. 

Two independent investigations of the ambient level 
of turbulence in the empty test section (with no test plate 
installed) were conducted. The normalized energy density 
spectra for these investigations are presented in Pig. 0-15. 

Both spectra were taken at the tunnel location corresponding 
to x/t = 3.25 and y/t = 0. For each investigation the value 
of u measured by the total and static head tubes was found to 
be 9.52 fps. Turbulence Intensities for these measurements 
were 0.0352 and 0.0366, considerably lower than the value 
of 0.068 found by Huval when operating at u = 10.2 fps. 

Also presented in Fig. C-15 are the theoretical 
functions suggested respectively by Dryden and von Karman: 

0(n) = -^1. and /(n) = ^2 

1 + B^n 2 (1 + B 2 2 n 2 )5/6 

While reasonable agreement is obtained with either function 

it appears that the experimental data are better represented 

by the Dryden expression in this case. The individual plate 

wake spectra are also well represented by this expression in 

the high frequency region. This is an indication of an 

approach to isotropic turbulence in this region, while anisotropy 

is Indicated in the low frequency plate spectra regions. 
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Pigs. C-1 through C-3 illustrate the downstream develop- 
ment of the turbulence spectra for the individual plates 

I 

along the centerline. Pigs. C-4 through C-6 present the 
longitudinal development for the individual plate spectra 
in line with the edge of the plate. Pigs. C-7 through C-10 
compare the centerline spectra for the three test plates 
at each of the four longitudinal stations, and Pigs. C-11 
through C-1 4 make a similar comparison for the spectra 
obtained at y/t = 0.5. Inspection of these Figures shows 
the following trends which are worthy of note: 

Between about 400 cps and 750 cps the slope of all 
the spectra is essentially -2 on the log-log plot. At 
frequencies between the double vortex shedding frequency 
and about 400 cps the slope is well represented by -5/3. 

Discrete energy density spikes associated with vortex 
shedding are very pronounced at x/t = 1.50 and at x/t = 

3.25 for both y/t = 0 and y/t = 0.5. Spikes associated 
with the double shedding frequency are pronounced only at 
the centerline. 

Considerable scatter exists in the spectra in the frequency 
range from 2.5 cps to about 20 cps. 

The discrete frequency energy appears to diffuse initially 
to the random components Immediately adjacent in frequency. 

These components in turn transfer part of their energy to 
their neighboring components. Thus, as the flow moves down- 
stream the discrete energy spikes are rapidly reduced by 
diffusion of their energy to the random components in a sort 
of "relaxation" process. 

As the flow moves downstream from x/t = 1.50, the 
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turbulence spectra indicate transition from very pronounced 
anisotropic turbulence superposed upon strong discrete 
components of vorticity toward a more isotropic condition 
of statistical equilibrium. Apparently the high frequency 
components achieve equilibrium early in the wake, since the 
portions of the spectra beyond 300 cps do not appear to 
undergo much change as the flow moves downstream. Hinze 
suggests a plausible reason for this early achievement of 
stability in the high frequency end of the spectrum. He 
points out that the flow within the smallest eddies is no 
longer turbulent, but viscous, and that the lower limit on 
turbulent eddy size is determined by dominance of molecular 
effects In general, the severity of velocity grad- 

ients increases with decreasing eddy size. The rate of 
viscous dissipation of energy associated with large velocity 
gradients then sets the lower limit on eddy size in accordance 
with the principles of conservation of energy. 

Numerical integration of the first four energy density 
spectra taken for Plate 0 resulted in very large energy defects 
compared to all the other spectra taken. Applying the 
correction suggested above as a possible means of adjusting 
for differences in measuring instruments (i.e., normalizing 
with respect to the Ballantine voltage recorded instead of 
the voltage reading of the thermocouple voltmeter) only 
slightly improved the situation. For these spectra A g still 
remained less than 50 % . 

Possible explanations for this large discrepancy are: 

(1) Incorrect wave analyzer calibration constant. This 
is discussed further in Appendix A. 
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(2) Abnormally large amounts of energy associated with 
frequencies above 750 cps. As shown in Fig. VIII, an extended 
spectrum analysis taken with the Type 1554-A analyzer indicates 
about 10$ of the spectral energy above 750 cps. However, the 
energy above 750 cps would have to be considerably larger than 
this to account for the discrepancy. 

(3) Presence of an air bubble in the probe tip. If 
the natural response of the gage-water-column-probe-tip 
system is located at a frequency just above the range of the 
wave analyzer, a large amount of energy will exist at this 
point and will not be measured by the analyzer. It will, how- 
ever, be measured by the rms voltmeter. This is considered 
unlikely, however, since two separate probe tip filling 
operations were involved. 

Fig. VIII indicates that the relative magnitude of the 
four Plate C spectra in question are probably correct, since 
the pattern is essentially duplicated by the spectrum obtained 
with the Type 1554-A analyzer. 

In order to determine whether the probe tip itself was 
subject to vibrations in the highly fluctuating flow field, 
the probe tip opening was plugged with candle wax. The 
turbulence probe was then located at x/t = 3.25 and y/t =0.5 
in the wake of Plate C. With a flow of U 0 = 9»55 fps, a 
turbulence intensity of 0.0019 was obtained. With the EinTel 
amplifier gain increased to 100, the entire spectrum was scanned; 
response was negligible at all frequencies. With the flow secured 
a voltage was still measured with the thermocouple voltmeter; 
this background signal, referred to a flow rate of 9.55 fps, 
corresponded to an intensity of turbulence of 0.0006. 
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A graphical presentation of the decay of discrete vortex 
energy is given below in Pig. IX: 





FIGURE IX Decay of Vortex Energy Downstream. 
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VIII. CONCLUSIONS 

A study of the graphical presentation of results Included 

as Appendices B and C Indicates the following trends for these 

4 

test plates at a plate-thickness Reynolds number of 1.9 x 10 : 

(1) Sharp geometry transitions (well defined flow 
separation points) in the trailing edge create a higher level 
of turbulence intensity in the early wake than do gradual 
transitions. Turbulence Intensity diffuses rapidly out of 
the wake at a rate proportional to the initial magnitude of 
the intensity and gradually approaches a decay law propor- 
tional approximately to (x/t)~£. 

(2) Strong wake underpressures exist in the early wake 
behind sharp geometry transition trailing edges. 

(3) The energy density spectrum in the early wake 

is anisotropic, particularly close to the trailing edge where 
longitudinal and transverse pressure gradients are greatest, 
and particularly for low frequencies. 

(4) The high frequency end of the energy density 
spectrum in the early wake is seen to approach isotropic 
conditions. This is also true of the ambient tunnel spectrum, 
E(n) is seen to vary essentially as n“ 2 in this region. 

(5) The discrete energy represented by the vortex 
shedding frequency is seen to decay downstream in a manner 
similar to the decay of turbulence intensity. The rounded 
trailing edge geometry favors the formation of vortices of 
higher energy content. 

(6) The magnitudes of the turbulence intensity found 
in the early wake close to the trailing edge are in excess 
of the order of magnitude originally assumed in the 



development of the expression for turbulence gage response. 
These magnitudes must, therefore, be viewed In a qualitative 
sense only. 

(7) A significant percentage (as much as 50/0 of the 
total kinetic energy of a mixed periodic-turbulent wake Is 
associated with discrete vortex shedding frequencies. 

(8) Energy associated with vortex shedding frequencies 

Increases Initially In the case of the plate with a 90° 

reentrant angle trailing edge (Plate A), and then decays 

(13) 

approximately linearly downstream. Huval found similar 
behavior in the early wake behind a 60° Included angle trail- 
ing edge geometry. 

(9) Taking the amount of energy associated with the 
discrete shedding frequencies in the immediate vicinity of 
the trailing edge as a measure of the energy available to 
excite plate vibrations, extrapolation of the trends indicated 
in Pig. IX indicate that Plate B would tend to vibrate most 

and Plate A would vibrate least. This is in agreement with 

' (14) 

results reported by Ippen et al 

While centerline vortex energy close to the trailing edge 

is greatest for Plate A, wake centerline forces exerted on the 

trailing edge will act equally on either side of the reentrant 

angle; response to such forces would thus inherently be damped 
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IX. RECOMMENDATIONS 

The following recommendations are made for the conduct 
of further investigation of the early wake behind flat plate 
structures : 

(1) The predictability of vortex shedding frequency and 
the behavior of plate-thickness and vortex-spacing Strouhal 
number should be investigated as a function of Reynolds 
number. Vortex spacing should be Investigated by traversing 
the early wake with the turbulence transducer and obtaining 
transverse distributions of discrete spectral energy. 

(2) The behavior of profile drag coefficient with var- 
iations In trailing edge geometry should be the subject of 

a separate investigation. Techniques and instrumentation 
in this investigation were Inadequate for any conclusive 
understanding of such behavior. 

(3) Particular care should be taken to duplicate both 
flow velocity and tunnel test section pressures for absolute 
comparison of successive runs. 

(4) Tunnel wall spindle cavities should be plugged with 
plates drilled to accept plate mounting spindles in order to 
eliminate the possibilities of Helmholtz or stub tube reson- 
ance effects. (See Appendix D. ) Any other cavities in the 
tunnel walls in the vicinity of the test plate or transducer 
should be similarly plugged or faired. 

(5) Frequent periodic calibration of measurement equip- 
ment and components!, particularly when they are battery 
powered, is mandatory when accurate data is to be obtained 
over a period of time. 
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(6) Accurate analysis of energy density spectra requires 
accurate calibration of analysis equipment. True rms indication 
of both total energy level and energy density as a function 

of frequency is desirable. The same instrument or accurately 
cross-calibrated instruments should be used to measure these 
quantities. 

(7) Consideration should be given to utilizing frequency 
modulated tape recording equipment to permit determination 

of autocorrelation functions for the energy density spectra. 
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APPENDIX A 

CALIBRATION PROCEDURE 
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. 1. Turbulence Pressure Transducer . 

a. General . 

The turbulence gage used in this investigation was the 

r n i 

same as that employed by Huval 1 JJ . A quasi-static calibration 
procedure similar to that used by Huval was used. However, 
the procedure was simplified as a result of repairs to the 
transducer made by Hydrodynamics Laboratory personnel which 
resulted in a greatly increased gage voltage decay time 
constant. 

Huval found the turbulence gage sensitivity to be 63 

millivolts per psi, a value substantially lower than the 

theoretically predicted sensitivity of 100 millivolts per 

T 13 1 

psi. Furthermore, his calibration 1 - ~ >J showed the gage time 
constant to be in the vicinity of 0.6 second, considerably 
lower than that to be expected. Subsequent to Huval's in- 
vestigation the sensitivity of the gage decreased markedly. 
Hydrodynamics Laboratory personnel then disassembled the 
gage and found a small amount of foreign matter on the 
ceramic face. The ceramic was carefully cleaned and re- 
assembled. Tests performed prior to the present investi- 
gation disclosed that the sensitivity of the gage was in- 
creased and the gage time constant was greatly improved. 

b. Description of Calibration Procedure . 

Fig. A-l shows schematically the turbulence gage cali- 
bration equipment used in this investigation. 
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FIGURE A-1 Gage Calibration Schematic 

The turbulence gage with probe tip installed was in- 
serted in the tube connected to the side of the cylindrical 
pressure tank as shown. A differential pressure was built 
up in the mercury manometer by means of the hand-operated 
pressure bulb. The clamp in the manometer pressure tube was 
set to give sufficient damping constriction to reduce pressure 
oscillations. The solenoid valve at the top of the pressure 
chamber was opened by means of a trigger circuit which simul- 
teneously triggered the horizontal sweep on the oscilloscope. 
An oscilloscope camera was used to record the indicated gage 
response voltage. See Fig. A-2. 

A simpler method of calibration, eliminating the use of 
the oscilloscope camera, was employed for the majority of the 




?ISUS 3 1-2 ixa-ple of Transducer- 
Preamplifier Gall brat ion Pro t oz rant. . 
(Photograph is reversed: sveep ori- 
gin is at the right and roves to 
the left. ) 
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turbulence gage sensitivity calibrations in this investigation. 

In this alternative method, the oscilloscope sweep rate was 
increased and the horizontal sweep selector set to "recurring 
sweep." The vertical deflection of the horizontal sweep at 
time of pressure release was measured directly on the face of 
the oscilloscope, manually recording on a strip of paper. 

One series of photographic calibration runs was made immediate- 
ly after a series of runs using the simplified procedure. 
Reproducibility between the two methods was excellent, both 
yielding the same value for gage sensitivity. 

c. Development of the Expression for Turbulence Gage 
Sensitivity and flme Constants . 

The following definitions will be used in this develop- 
ment ; 

K^g = turbulence gage sensitivity 

Ktp = transducer-preamplifier system sensitivity 

Gp = preamplifier gain 

T v = voltage time constant 

Tp = pressure time constant 

C v = 1 /T v 

c p = 'Ap 

p = Instantaneous pressure at transducer 

p 0 = initial calibration pressure 

V = instantaneous voltage at oscilloscope 
As Indicated by Huval^^, the sudden release of 
pressure from the tank, assuming a suitably damped system, 
may be represented by the following expressions; 
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p = PqC - ^ , for t >0 

P = P Q » for t <0 (A-l) 

The incremental change in voltage generated by the tur- 
bulence gage is the sum of two components: 

AV]_ = -Kf-pAp (pressure-generated component) 

AV 2 = -C v VAt (voltage decay component) 

Hence: AV = AV 2 + AV 2 = -K tp Ap - C y VAt (A-2) 

Rewriting Eq. (A-2) as: 

= ' C v VAt 

and taking the limit as At -► 0, we obtain: 

M = " K tp at ' C v V < A -3> 

Substituting p = Poe’^P^ into Eq. (A-3) we obtain 

3F = K tp p o°p e ' Cpfc * °v v < A -W 

Defining V Q = K tp p Q , we see that: 

HF = V o C p e ' Cpt - c v v (A-S) 

Multiplying both sides of Eq. (A-£) by e^ v ^: 

e ° vt HF + c vVe Cvt = V 0 C p e( C v” C p) t (A-6) 

Then: 

4- (v e Cvt ) = V 0 C p e' C P fc (A-7) 

Integrating: 

V e c vt = VoCp e (Cv-C p )t + 



(A-8) 
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where C 2 is the constant of integration. From the initial 
condition that V = 0 at t = 0, we find that: 



Then: 




(A-9) 



(A-10) 



In analysis of the oscilloscope calibration photo- 
graphs, it is convenient to express V in terms of V m , the 
maximum voltage occurring at time t m . (See Fig. A-2). 
From: 



dV 



t = t, 



= 0 



m 



we obtain: 



Also: 



e" C P t m = e “ Cvtm 
C P 



(A-ll) 



v m = !2gp ( 

c v" c p \ 



e" C P t m _e" C v t m 



( A-12 ) 



Substituting Eq. (A-ll) into Eq. (A-12) and simplifying: 

v m = Voe-^m = KtpPoe-W^v (A . 13) 

Solving for K^ p in terms of V m » we obtain: 

_ V m W T v (A-]4) 

K tp - e 

Before going on to the development of expressions for 

the time constants T v and Tp, it should be noted from the 

foregoing expression that if t m <<^ T v , then: 

K t n * ^SL 

tp Po 



(A-15) 
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This yields a means of obtaining an adequate approx- 
imation for without accurate knowledge of either T y 
or T p . 

Substituting Eq. (A-13) into Eq. (A-10), the following 
expression for V may be found: 



V = 






e 




-e 





(A-16) 



If it can be assumed that T p « T v and t m «T v , the 
second term within brackets in Eq. (A-16) becomes negligible 
with respect to the first term for t » t m . 

Then: 

- \: t|n 

(A-17) 



V m e 



t , 7 t m 



} for t » t 



m 



Hence: 

T v = ln (Vm/v) ; evaluated at t^> t m (A- 18 ) 



To obtain a means of evaluating T p , evaluating dV/dt 
from Eq. (A-16) above at t = 0 results in: 



and: 




= 0 



_ v m 





V m . 

[dV/dt J 



; evaluated at t = 0 (A- 19 ) 



This allows evaluation of T p from knowledge of the 
initial slope of the V versus t plot. 

d. Results of Pressure Transducer Calibration. 



A series of seven photographic calibration runs were 
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made using various values of initial pressure p 0 and using an 
oscilloscope horizontal sweep rate, calibrated, of $00 milli- 
seconds per inch. A DuMont Oscillograph-Record camera Type 
297, using an f-stop of f/5.6 with a 5-second exposure, the 
sweep being triggered 2 seconds after the shutter was opened. 
Polapan 3000 Type 14-7 film was used. The time to reach maximum 
voltage, t m , was about 30 milliseconds on each run, as com- 
pared with a value of about $0 milliseconds reported by Huval. 

(1) Pressure Decay Time Constant, Tp . 

Eq. (A-19) was used to determine Tp. Values ranged 
from 11 to 18 milliseconds, with five of the seven values 
being in the range of 12.5 to 1J4. milliseconds. An average 
value of Tp = 13*7 milliseconds was obtained. 

(2) Voltage Decay Time Constant, T v . 

A plot of ln(Vm/v) versus (t-t m ) was used to 
evaluate T p . See Eq, (A-18). A range of values from 2,i4. 
seconds to 11.1+ seconds was obtained, with five of the seven 
runs in the range I+.I4. to 9.0 seconds. 

Mubh of this scatter was due to inability to obtain an 
absolutely stable reference voltage on the oscilloscope prior 
to triggering the pressure release. Also, extreme sensitivity 
of the gage and preamplifier to small unavoidable disturbances 
during the 2-second period of the horizontal sweep may have 
contributed to the scatter. 

The two extreme values were disregarded, resulting in 
an average from the remaining five runs of Tp = 6.I4. seconds. 

However, even using the smallest value of T y and the 
largest value of T p , it is seen thats 
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This shows that the approximation of Eq. (A-l£) is 
valid for determining K tp and justifies using the simplified 
procedure described in Section A.l.b above. 



Seven separate series of calibration runs were performed 
during the course of this investigation between 22 February 
and 11 April 1961. The dates of the runs are presented in 
the legend of Fig. A-3. 

During the observations made on 6 April, the KinTel 
amplifier was introduced into the calibration circuit to 
provide higher output voltage levels. Amplifier gain setting 
was 30. It was felt that this procedure would serve as a 
further check on the sensitivity of the entire measuring 
system. Due to the quasi-static nature of the calibration 
procedure, the DC filter between the preamplifier and the 
amplifier was bypassed. The average value of K^. determined 
for these runs was about 20°/o lower than the February cali- 
bration. 

As a check on possible time-deterioration of sensitivity 
the February procedure was duplicated on 11 April; these runs 
produced a value of about l±°/o lower than the original 
February value. It was concluded that the amplifier musu 
have attenuated the quasi-steady gage response. 



(3) Turbulence Gage Sensitivity . 

Turbulence gage sensitivity was determined from 



Eq. (A-15) and the relationship: 
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FIGURE A-3 Transducer-Preamplifier Sensitivity Calibration 



PRESSURE DIPPEElENTIAL , Ah , FEET OP WATER 
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The amplifier had already been calibrated using an AC 
signal source and found to have excellent gain response. The 
amplifier was now calibrated using a DC dry cell through a 
potentiometer as a signal source. With gain setting of 30, 
actual gain was measured to be equal to about 27. When the 
ratio of gain factors was applied to the gage calibration 
constants obtained with the amplifier in the circuit, the 
same b°/o variation from the initial calibration sensitivity 
was obtained. It is considered probable that this variation 
is due to data scatter. 

Finally, an average value of = 16.0 mv/ft. of water 
= 36.9 mv/psi was obtained. This corresponds to a value of 
= i+O mv/ft. of water = 92.1]- mv/psi, as compared to the 
theoretical prediction of = 100 mv/psi reported by 
Perkins and Eagleson tJ . 

2. Wave Analyzer . 

Both wave analyzers used in the investigation were 
calibrated using a "pure tone" sinusoidal input signal from 
the Hewlett-Packard audio oscillator. Rms input voltages 
were measured with the Ballantine VTVM. All calibrations were 
performed with the vibration analyzer selectivity in the 
"sharp" position and with vibration analyzer sensitivity set 
at 10. (All data were taken with these settings also.) 

a. General . 

As Indicated by Perkins and Eagleson^^, the General 
Radio Type 762-B Vibration Analyzer output voltmeter reading 
squared is proportional to the total energy level of the 
frequency components of the Input signal within the narrow 
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band to which the analyzer is tuned. 



b. Wave Analyzer Response . 



^bw 



t An 

“c+T 



~ An 
n c“T 



E(n)dn 



(A-20) 



For an input signal with a continuous spectrum of 
frequencies, if E(n) may be considered constant over the 
narrow bandwidth of the wave analyzer, (A-20) becomes: 

E bw = E(n c )An = C]_ E(n c )n c (A-2l) 

From the statement made in paragraph (a) above, there 
results : 

E bw = C 2 Rd 2 (A-22) 

Therefore, combining (A-2l) and (A-21) : 

E(n c ) = (A-23) 

°v n c 

C 1 

where C y = -^ 

c. Wave Analyzer Calibration Procedure . 

The purpose of the wave analyzer calibration is to de- 
termine the constant C v in (A-23) • 

One means of obtaining this constant is to set the wave 
analyzer at a given frequency n^, then to apply an input 
sinusoidal signal from a sine-wave generator. The input 
signal frequency is increased in small Increments from a 
value (well below n^) where negligible wave analyzer response 
is noted, to a value well above n^ where again negligible 
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analyzer response results. The input rms voltage e^, as 
measured by an A.C. rms voltmeter, should be kept constant 
throughout this procedure at a value predetermined to give 
nearly full scale deflection at maximum wave analyzer response, 
e^ is recorded and constantly monitored. At each input fre- 
quency n, n and Rd are recorded. Then C v may be determined 
from the relationship; 

~ca 

.2 



c = -A 

v 

e i 



Rd~ 
r i. 



dn 



(A-2l*) 



The foregoing procedure should be repeated for several 

values of n over the range of frequencies which are of 
c 

interest. By this means, any frequency dependency of C v may 
be detected. 

An alternative procedure (introducing a small error) 
is to set the input frequency at a constant value n-^, then 
shift the wave analyzer frequency in small increments from 
the point of negligible response on the lower side of n^ to 
the point of negligible response on the upper side of n-, . 
Then C y may be approximated from the relationship; 



pCD 



A —o ^ 

c v ~ “T 
e i 



Rd 2 



n, 



dn. 



e 



(A~2£) 



The latter method was used throughout this investigation, 
due to instability of the sine wave generator output voltage 
introduced by shifting frequencies. Numerical integration 
was used to approximate the value of the integral. 
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d. Results of Wave Analyzer Calibrations . 

It was noted during this investigation that the wave 
analyzer calibration constant was dependent upon plate vol- 
tage (supplied by two i|£ volt batteries). 

This was true even when a check of plate battery voltages 
showed the plate voltage to be within the satisfactory 
operating range. Therefore, frequent wave analyzer calibration 
is necessary in order to obtain accurate results. Particular- 
ly, the wave analyzer must be recalibrated whenever the 
plate supply batteries are changed. 

It is possible that large discrepancies in the areas of 
four energy spectra for Plate C may be attributed to failure 
to calibrate the wave analyzer immediately after a change of 
plate supply batteries. 

On March 28, 1961, wave analyzer plate supply batteries 
were replaced. The turbulence spectrum measurements were 
commenced and all eight energy spectra for Plate C were 
taken between March 28 and March 31* when the wave analyzer 
was calibrated. The wave analyzer constant was found to be 
in good agreement with the two previous values determined. 
However, it was later discovered that the integrals of the 
energy density spectra were very low for the four spectra 
measured on March 28 and March 29 compared to those taken on 
March 30, and. In fact, compared to all the spectra taken for 
the other two plates. At present, lacking a better reason 
for this discrepancy, it is conceivable that the wave analyzer 
constant was very low Immediately after the battery replacement 
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for some unexplainable reason. In evaluation of this possi- 
bility, it must be noted that a very low value compared to 
any previously obtained would have been necessary to explain 
the discrepancy, and the constant would have increased with 
analyzer operating time, contrary to the trend previously 
experienced. 

The following values of C v were obtained during this 
investigation. Only average values are listed, since each 
calibration showed little frequency dependence of the constant. 



3. Amplifier and Preamplifier . 

The electrometer preamplifier and the KinTel amplifier 
were both calibrated using a sinusoidal input signal from 
a Hewlett-Packard Model 2000 audio oscillator. 

Calibration of these components was conducted before 
and after data takings no time variation of response was 
observed. The Initial calibration was performed for both 
instruments by measuring both input and output signals with 
the Ballantine VTVM. This necessitated removing the volt- 
meter from the input circuit to measure output voltages from 
the component being calibrated. While this procedure is not 
considered sound practice, the results of the initial cali- 
bration were corroborated by subsequent calibrations. The 
low frequency limit of the audio oscillator used as a signal 
generator is 16 cps. Hence, the equipment could not be 



PERIOD 




3-6-61 to 3-10-61 
3-10-61 to 3-28-61 
3-28-61 to 3-31-61 



5.12 x 105 
£.£1 x 105 

5.31 x io£ 
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calibrated for frequencies below this limit. The initial 
calibrations were conducted from 16 cps in harmonic increments 
of 60 cps to an upper limit of 720 cps with continuous oscil- 
loscope monitoring. Gain factors were computed for the pre- 

♦ 

amplifier, and for the amplifier at gain settings of 20, 30, 
£0, 70, and 100. 

The final calibrations were conducted using the Ballan- 
tine VTVM for input measurements and the thermocouple volt- 
meter for output measurements. Simultaneous measurements were 
taken without disturbing the calibration circuitry. The two 
meters were previously calibrated against each other with a 
sinusoidal input signal over the entire frequency range. 
Finally, the preamplifier and amplifier were calibrated to- 
gether as a system from 16 cps to £000 cps with amplifier 
gain settings of 20 and 30 (i.e., the gain settings actually 
used during data taking). 

The results of these calibrations are presented in Fig. 
A-4 as a plot of relative attenuation or loss in gain in db 
versus frequency. 

The measured gain of the preamplifier is seen to exhibit 
some attenuation above 350 epsj the response of the amplifier 
is essentially flat up to 750 cps, which is the upper limit 
of the Type 762-B Vibration Analyzer. The preamplifier gain 
was taken as constant over the operating range of the in- 
vestigation at a value G - Q.ij.Q, which was the arithmetic 

s' 

mean value. The amplifier gain was G_ = 20 with a gain 

a 

setting of 20 and G a = 30 with a gain setting of 30. 
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FIGURE A-4 Component and System Gain Response 



The electrometer preamplifier was calibrated and 
operated throughout with a screen grid bias of +2.0 volts 
A circuit diagram of the electrometer preamplifier is pre 



sented in Pig. A-5. 
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FIGURE A-5 Electrometer Preamplifier Circuit Diagram 
(From Huval, "The Early Wake of a Fixed, Flat Plate.") 
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FIGURE B-1 Wake Structure Summary Plot, Plate A 
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FIGURE B-3 Wake Structure Summary Plot, Plate 0 
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FIGURE B-4 Transverse Velocity Profile, Plate A 
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FIGURE B-5 Transverse Velocity Profile, Plate B 
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FIGURE B-6 Transverse Velocity Profile, Plate 0 
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FIGURE B-7 Comparison of Velocity Profiles at x/t = 1„50 
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FIGURE B-9 Comparison of Velocity Profiles at x/t =6.50 
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FIGURE B- 10 Comparison of Velocity Profiles at x/t = 9„75 
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FIGURE B-12 Turbulence Intensity Distribution at x/t = 
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FIGURE B-15 Pressure Coefficient Distribution at x/t = 
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FIGURE B-20 
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FIGURE B-25 Wake Symmetry Plot, Plate 0 
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APPENDIX C 

NORMALIZED ENERGY DENSITY SPECTRA 
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FIGURE 0-1 Energy Density Spectra at y/t = 0, Plate A 
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FIGURE C-2 Energy Density Spectra at y/t = 0, Plate B 
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FIGURE 0-3 Energy Density Spectra at y/t = 0, Plate 0 
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FIGURE 0-4 Energy Density Spectra at y/t =0.5, Plate A 
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FIGURE 0-5 Energy Density Spectra at y/t =0.5, Plate B 
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FIGURE 0-6 Energy Density Spectra at y/t = o.5, Plate C 
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FIGURE 0-7 Energy Density Spectra at y/t = 0, x/t = 1 „ 50 
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FIGURE 0-8 Energy Density Spectra at y/t = 0, x/t = 3o25 
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FIGURE C-9 Energy Density Spectra at y/t = 0, x/t =6.50 
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FIGURE 0-10 Energy Density Spectra at y/t = 0, x/t = 9.75 
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FIGURE 0-11 Energy Density Spectra, y/t = 0.5, x/t = 1.50 
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FIGURE C-12 Energy Density Spectra, y/t = 0.5, x/t = 3.25 
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FIGURE 0-13 Energy Density Spectra, y/t = 0.5, x/t = 6.50 
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FIGURE 0-1A Energy Density Spectra, y/t = 0.5, x/t = 9.75 
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FIGURE C-1 5 Energy Density Spectrum, Empty Test Section 
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APPENDIX D 

HELMHOLTZ RESONANCE EFFECTS 



1. Plate Spindle Cavities . 

During measurements incident to determination of the 
characteristics of the ambient tunnel turbulence spectrum 
a possible Helmholtz resonance effect was observed. With 
the desired flow established, the amplified turbulence in- 
tensity voltage was observed to be lj.0 mv rms. Monitoring 
the entire available spectrum with the vibration analyzer 
showed no significant peaks. The test plate had been re- 
moved from the tunnel test section; no plugs had been placed 
in either of the spindle cavities. 

After several minutes the turbulence voltage level began 
to rise, reaching approximately 90 mv rms, and a discrete 
quantity of energy was observed with the vibration analyzer 
at about lj.0 cps. The remainder of the spectrum indicated an 
essentially constant level, down about 6 db from the peak 
reading. The rms voltage level gradually increased and the 
peak frequency was gradually noted to decrease, reaching as 
low as 33 cps. 

At this point a large accumulation of air bubbles was 
noted in the upper spindle cavity under the top cover plate. 
The vent was opened, and as the bubbles vented off, the rms 
voltage level decreased to about i;0 mv and the discrete peak 
disappeared. 

This phenomenon indicated the possibility of Helmholtz 
resonance to the authors. Clearly only the air bubbles were 
causing the low frequency energy, but even a water-filled 
cavity might be expected to resonate and contribute energy to 
the flow. Even with the test plate installed and the mounting 



spindles passing through the upper and lower tunnel walls, 
cavities of complicated geometry exist. See Fig. D-l below. 



1 24 




The fundamental resonant frequency of a Helmholtz type 

rru 

cavity can be estimated using the electrical analogy of 
a simple tank circuit. For a cavity whose dimensions are 
known: 

S = surface area of the cavity opening, 

V = volume of the cavity, 

l e = effective length of the "neck" of the cavity, 
approximated by actual length + 0.8y^, and 
c = sonic velocity in the medium, taken as k-79b 
fps in fresh water, uncorrected for rigidity 
of cavity "neck. " 



Then the Hplmholtz resonant frequency, njj, is approximated by: 



and in this case from actual measurements: 

n H = ^Zw * 1 y = 26 ^° cps 

It is well known that the theoretically derived 
frequency is an upper limit and that in general any number 
of effects can reduce the derived frequency. Furthermore, 

rru 

Morse points out lJJJ that the lumped circuit approximation 
is less valid for acoustical systems than it is for electrical 
circuits because of the shortness of the wave lengths in- 
volved. He notes that Eq. (D-l) is of questionable validity 
for frequencies above 1,000 cps, for when a system is longer 
than a half wave length it begins to behave like a trans- 
mission line and the wave equations become involved. 

2. The Turbulence Probe Tip as a Helmholtz Resonator . 

The above discussion led naturally to a consideration of 
the turbulence probe tip as a Helmholtz resonator. From the 
interior dimensions of the probe tip as cited in Reference (13) 
and using Eq. (D-l), it was found that for the turbulence 
probe tips 

n H ~ 7060 cps 

It is considered probable that signals above this fre- 
quency will be attenuated at the rate of approximately 20 db 
per decade. 
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APPENDIX E 

REPRESENTATIVE SAMPLES OP DATA 
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SAMPLE WAKE SURVEY DATA 



PLATE 0. x/t = 9.75 DATE: 3-28-61 

AMPLIFIER GAIN: 50 



CONTRACTION MANOMETER: Ah M = 21.30, 


Ah HR =1.27 


VALVE INDICATORS - FOOT VALVE: 


1 1 9 


PUMP DISCHARGES 


22.5 


"CIRCUIT STOP: 





TURBULENCE VOLTAGE (RMS) TOTAL HEAD STATIC HEAD 



y/t 


Thermocouple 

mv 


Ballantine 

mv 


Pressure 

Side 

"Hg 


Atmos. 

Side 

"Hg 


Pressure 

Side 

"Hg 


Atmo s . 
Side 
"Hg 


+1.5 


1 10 


93 


20.73 


25.86 


20.55 


24.47 


+1 .4 


115 


102 


20.74 


25.85 


20.56 


24.47 


+1.3 


120 


105 


20.75 


25.84 


20.56 


24.46 


+1.2 


125 


1 10 


20.76 


25.84 


20.56 


24.46 


+1 . 1 


135 


122 


20.78 


25.82 


20.57 


24.45 


+ 1 .0 


139 


127 


20.80 


25.77 


20.57 


24.45 


+0.9 


145 


130 


20.82 


25.75 


20.57 


24.45 


+0 . 8 


151 


132 


20.83 


25.74 


20.57 


24.45 


+0.7 


153 


136 


20.84 


25.73 


20.58 


24.44 


+0.6 


156 


138 


20.86 


25.72 


20.58 


24.44 


+0.5 


157 


138 


20.87 


25.71 


20.58 


24.44 


+0 .4 


155 


142 


20.88 


25.70 


20.58 


24.44 


+0.3 


155 


1 4 1 


20.90 


25.68 


20.59 


24.43 


+0.2 


157 


143 


20.90 


25.68 


20.59 


24.43 


+0.1 


154 


142 


20.90 


25.68 


20.59 


24.43 


0 


150 


138 


20.90 


25.68 


20.60 


24.42 


-0. 1 


155 


143 


20.90 


25.68 


20.59 


24.43 


-0.2 


157 


145 


20.88 


25.70 


20.59 


24.43 


-0.3 


154 


1 40 


20.86 


25.73 


20.59 


24.43 


-0.4 


149 


142 


20.86 


25.73 


20.59 


24.43 


-0.5 


149 


139 


20.82 


25.76 


20.57 


24.45 


-0.6 


150 


138 


20.80 


25.78 


20.56 


24.46 


-0.7 


147 


132 


20.80 


25.78 


20.56 


24.46 


-0.8 


142 


130 


20.77 


25.81 


20.56 


24.47 


-0.9 


135 


123 


20.76 


25.83 


20.56 


24.47 


-1 .0 


133 


123 


20.75 


25.84 


20.56 


24.47 


-1 . 1 


122 


110 


20.74 


25.85 


20.55 


24.47 


-1.2 


1 17 


103 


20.74 


25.85 


20.55 


24.47 


-1.3 


115 


95 


20.73 


25.86 


20.55 


24.47 


-1 .4 


100 


87 


20.72 


25.86 


20.55 


24.47 


-1.5 


97 


82 


20.72 


25.87 


20.56 


24.48 



SAMPLE ENERGY SPECTRUM DATA 



PLATE C. 


x/t = 3.25. 


v/t = 0.5 






DATE: 3-29-61 


AMPLIFIES 


: GAIN: 20 










AVERAGE RMS TURBULENCE VOLTAGE 


- THERMOCOUPLE: 131 


mv 








BALLANTINE: 


115 


mv 


ANALYZER 


SENSITIVITY 


: 10 








ANALYZER 


SELECTIVITY 


: SHARP 








ANALYZER 


CALIBRATION 


CONSTANT: 


5.31 x 105 volts” 1 




CONTRACTION MANOMETER: Ahu = 21 


.35, Ah P „ = 1 


.?o 




Ail total = 


: 25.62 - 20 


.95 ^static = 2 4.36 


- 20. 


66 


FREQUENCY ANALYZER 


FREQUENCY 


ANALYZER FREQUENCY ANALYZER 


(cps) 


READING 


(cps) 


READING 


(cps) 


READING 


n 


Rd 


n 


Rd 


n 


Rd 


2.5 


7.5 


70 


22 


230 


28 


3 


9.4 


75 


25 


240 


26 


4 


8.9 


80 


28 


250 


24 


5 


9.5 


35 


31 


260 


25 


6 


1 1 .0 


90 


38 


270 


24 


7 


9.9 


95 


60 


280 


27 


8 


11.7 


98 


75 


290 


26 


9 


12 . 1 


100 


92 


300 


27 


10 


12.5 


102 


95 


320 


25 


1 1 


12.6 


104 


88 


340 


25 


12 


1 2.6 


106 


75 


360 


25 


13 


12.4 


108 


70 


380 


25 


14 


12.4 


110 


60 


400 


25 


15 


12.5 


115 


45 


420 


24 


17 


13.4 


120 


39 


440 


23 


20 


1 1.8 


130 


30 


460 


22 


22 


12.2 


140 


28 


480 


22 


25 


1 2.8 


150 


26 


500 


22 


28 


14.0 


160 


27 


520 


22 


30 


14.8 


170 


27 


540 


22 


35 


14.8 


180 


26 


560 


22 


40 


17.0 


190 


28 


580 


21 


45 


17.8 


195 


30 


600 


21 


50 


19.6 


200 


33 


650 


20 


55 


20 


205 


32 


700 


20 


60 


21 


210 


31 


750 


20 


65 


20 


220 


30 



















NOTE: Analyzer readings from n = 2.5 through n = 50 are time 

averages of a number of individual readings taken at fixed 
time increments in an attempt to smooth out data scatter. 
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